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£/>X>fk1-- 1^7*5.  ABSTRACT 


The  Transonic  Turbine  Test  Rig  installed  in  the  Turbo-Propulsion 
Laboratory  at  the  Naval  Postgraduate  School  was  designed  such  that 
losses  in  nozzles  and  rotor  blading  and  overall  performance  charac- 
teristics can  be  determined  for  operation  at  several  combinations  of 
blade  configurations,  pressure  ratios  and  speeds.   This  study  relates 
the  empirical  data  for  stator  loss  coefficients  to  those  determined 
for  the  converging-nozzle  stator  presently  installed,  and  discusses 
reasons  for  some  of  the  discrepancies.   Also  discussed  are  results  of 
the  flow  nozzle  calibration,  locked  rotor  data,  and  temperature  and 
inlet  swirl  effects.   A  method  is  proposed  to  reduce  these  effects 
for  future  work  on  the  Turbine  Test  Rig. 
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I.   INTRODUCTION 

Much  of  the  design  and  analysis  of  modern  turbines  is  based  on 
the  experience  of  the  designer  and  his  use  of  empirical  data  available 
in  the  literature.   As  a  re=iJ.t,  many  investigations  of  test  turbines 
are  being  carried  out  in  an  effort  to  more  fully  explain  and  mathe- 
matically model  the  flow  phenomena  in  stationary  and  rotating  channels, 
Present  methods  for  design  of  axial  flow  machines  base  calculations 
on  a  reference  diameter,  usually  taken  as  the  mean  of  the  hub  and  tip. 
A  major  simplification  is  then  made  by  assuming  that  at  any  point  in 
a  cross-section,  such  parameters  as  total  pressure  and  temperature  and 
axial  velocity  have  the  same  values  as  at  the  reference  diameter.   The 
overall  characteristics  of  a  stage  are  then  determined,  but  depend  on 
how  accurately  flow  discharge  angles  and  loss  coefficients  are  known. 
Thus,  it  is  extremely  important  to  be  able  to  predict  these  values  as 
functions  of  blade  geometry,  Mach  number  and  certain  other  flow  para- 
meters.  Such  a  method  as  just  described  is  conventionally  known  as 
a  "One-Dimensional  Analysis". 

The  Transonic  Turbine  Test  Rig  (TTTR)  located  at  the  Turbopropul- 
sion  Laboratory  is  a  particularly  desirable  test  machine  because  it 
permits  evaluation  of  performance  at  many  combinations  of  operating 
conditions  and  configurations.   One  of  the  most  important  features 
of  the  TTTR  is  that  it  enables  one  to  evaluate  flow  angles  and  loss 
coefficients  using  momentum  and  moment  of  momentum  equations,  without 
making  extensive  flow  surveys  which  cause  an  upset  of  the  true  flow 
conditions.   The  empirical  methods  available  for  prediction  of  flow 
angles  and  loss  coefficients  yield  results  varying  from  unreasonably 
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low  values  to  one  so  conservative  it  would  be  unlikely  the  machine 
would  operate.   Results  from  previous  runs  on  the  TTTR  showed  the 
latter  trend  in  loss  coefficients,  although  the  TTTR  had  a  reasonable 
overall  performance. 

This  thesis  investigates  the  methods  by  which  the  loss  coefficients 
are  determined  for  the  TTTR  and  relates  these  coefficients  to  those 
predicted  by  the  literature.   It  then  describes  the  investigation  car- 
ried out  to  explain  the  apparent  discrepancies,  and  examines  the  effect 
of  some  of  the  more  important  parameters  used  to  evaluate  the  loss  co- 
efficient. 

The  author  would  like  to  express  his  appreciation  for  the  compe- 
tent assistance  provided  by  Mr.  Jim  Hammer  and  members  of  the  Turbo- 
propulsion  Laboratory.   Special  thanks  are  extended  to  Dr.  M.  H.  Vavra 
for  his  guidance  and  advice  during  the  investigation. 
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II.   EQUIPMENT  MODIFICATIONS 

The  TTTR  instrumentation  has  been  described  extensively  by  Commons 
[Ref.  4]  and  Lenzini  [Ref.  6].   Discussion  here  will  be  limited  to 
recent  modifications  made  to  the  TTTR. 

Fig.  1  shows  the  piping  and  valve  location  in  the  Turbine  test 
cell.   Fig.  2  shows  a  cross-section  of  the  floating  stator  and  sup- 
port assembly,  while  Fig.  3  shows  a  more  detailed  cross-section  of 
the  stator  plenum  and  hub  area.   The  closure  plate  assembly  shown  in 
Fig.  3  was  instrumented  with  additional  strain  gauges  on  the  flex- 
ures, and  the  back  face  towards  the  rotor  had  three  static  pressure 
taps  installed  at  known  radii.   Results  of  these  pressure  readings 
and  closure  plate  torque  calibration  are  discussed  later. 

Fig.  4  shows  the  floating  stator  assembly  with  the  locations  of 
the  axial  force  and  torque  capsule  supports  and  Fig.  5  shows  a  more 
detailed  view  of  the  torque  capsule  support  assembly.   The  capsule 
support  assembly  was  moved  from  the  previous  location  between  the 
overhead  support  and  plenum  casing  to  its  position  at  the  upstream 
end  of  the  stator  in  an  attempt  to  reduce  both  the  temperature  effect 
on  the  capsule  and  the  errors  induced  by  differential  thermal  expan- 
sion. 

Fig.  6  shows  phe   stator  exit  and  shroud  area.   Two  additional 
pressure  taps,  P13  and  P14  were  installed,  with  the  hope  of  corre- 
lating the  stator  discharge  tip  pressure  to  the  shroud  pressure  at 
the  same  plane.   In  addition,  the  machine  was  completely  dismantled 
and  cleaned,  and  new  pressure  lines  were  installed  in  many  cases. 
The  pressure  lines  were  connected  to  a  mercury  manometer  board  in  a 
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logical  sequence,  extending  from  the  nozzle  to  the  plenum  assembly 
and  along  the  shroud  to  the  rotor  discharge,  so  that  the  pressure 
distribution  through  the  machine  could  be  readily  observed. 

Fig.  7  shows  the  converging  nozzles  presently  installed  in  the 
TTTR  and  Fig.  8  shows  the  circular-arc,  sharp  leading-edge  rotor  used, 
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III.   FLOW  RATE  DETERMINATION 


A.   NOZZLE  CALIBRATION 

Previous  nozzle  calibration  techniques  have  been  described  by 
Naviaux   and  Lenzini  [Ref.  6;.   It  was  shown  that  the  nozzle  coeffi- 
cient was  a  function  of  nozzle  Reynolds  number.   However,  former  re- 
sults for  nozzle  coefficients  were  determined  as  a  function  of  the 
Reynolds  number  based  on  the  conditions  of  the  ASME  standard  calibrat- 
ing orifice,  rather  than  on  the  flow  nozzle.   As  a  result,  it  was  de- 
cided to  repeat  the  nozzle  calibration  to  obtain  a  new  relationship, 
and  to  determine  if  any  deterioration  of  the  nozzle  had  taken  place. 

The  program  for  the  data  reduction  is  given  in  Appendix  D,  program 

2 
FLOCAL .   The  method  used  is  similar  to  that  given  by  Williams   ,  which 

was  found  to  be  more  suitable  for  computer  computations.   The  basic 

3 
equation  is  given  in  the  ASME  Power  Test  Code   and  values  of  the  con- 

4 
stants  are  taken  from  Stearns.  Tests  were  carried  out  at  nozzle  sup- 
ply pressures  of  26,  29,  34  and  39  psia  over  a  range  of  Reynolds 
numbers  from  2  x  10   to  12  x  10  .   Fig.  11  shows  the  results  of  these 


Naviaux,  J.  C,  "Transonic  Turbine  Test  Rig  Exhauster  System  Tests 
and  Tests  of  a  Reaction  Turbine,  "Naval  Postgraduate  School  Thesis, 
December,  1966. 

2 
Williams,  D.  D.,  'Determination  of  Performance  Parameters  of  a 

Dual-Discharge  Radial  Turbine,  "Naval  Postgraduate  School  Thesis, 

December,  1968. 

3 

Flow  Measurement,  Chapt.  4,  Part  5,  Supplement  to  ASME  Power  Test 

Codes,  ASME,  New  York,  N.  Y.,  1959,  p.  57ff. 

4 
Stearns,  R.  F.,  et_  al. ,  Flow  Measurement  With  Orifice  Meters. 

New  York:  D.  van  Norstrand  Company,  1951.  pg.  65-212. 
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tests,  with  nozzle  coefficient  plotted  as  a  function  of  Reynolds 
number.   The  curve  differs  from  that  given  by  Lenzini  [Ref.  6],  being 
higher  by  2  to  3  percent.   The  curve  also  shows  a  maximum  coefficient 
occuring  at  a  Reynolds  number  of  6  x  10   and  then  decreasing  as  Rey- 
nolds number  is  increased.   It  is  suspected  that  this  increasing  trend 
of  nozzle  coefficient  each  time  a  calibration  is  performed  indicates 
some  deterioration  of  the  nozzle. 

An  analytical  curve  for  the  nozzle  coefficient  was  obtained  by  us- 
ing a  least-squares  polynomial  curve  fit.   A  fourth-order  polynomial 
was  chosen  as  the  best  fit  to  the  data,  namely: 

C  =  9.5208828  x  10_1  +  4.4891640  x  10~7  Re  -7.5810743  x  10~13  Re  2 
n 

+  5.5023778  x  10~19  Re  3  -1.4571314  x  10~25  Re  4 

The  maximum  error  between  the  analytical  curve  and  the  actual  data 

5  5 

points  for  a  range  of  Reynolds  numbers  between  4  x  10   and  12  x  10   is 

0.2  percent. 

B.   TURBINE  FLOW  RATE 

The  program  to  obtain  the  turbine  flow  rate  is  given  in  subroutine 
FLORAT  of  program  TTTR,  and  employs  a  procedure  similar  to  that  de- 
scribed in  section  A.   The  labyrinth  leak  rate  uses  a  formula  deter- 
mined by  Lenzini  [Ref.  6]  which  is  a  function  of  plenum  pressure  ratio 
and  temperature.   The  labyrinth  leak  rate  was  found  to  be  only  2  per- 
cent of  the  total  flow  rate.   The  recalib ration  of  the  nozzle  and  im- 
proved flow  rate  determination  resulted  in  referred  turbine  flow  rates 
which  were  3  percent  higher  than  those  previously  used. 
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IV.   ANALYSIS  AND  DATA  REDUCTION 

A.   GENERAL 

The  basic  equations  for  a  one-dimensional  flow  analysis  of  the 
TTTR  are  well  known  and  have  been  described  in  detail  by  Commons 
[Ref.  4],  Lenzini  [Ref.  6]  and  Messegee  [Ref.  7].   The  reference  diam- 
eter referred  to  previously  is  twice  the  mean  radius: 


R  = 


R  ,  +  R  . 
hub     tii 


m        2 

The  data  for  the  TTTR  are  reduced  by  program  TTTR  given  in  Appendix  D, 
which  also  gives  the  results  of  Run  15.   Fig.  12  shows  the  thermody- 
namic process  for  an  axial  turbomachine  and  the  stator  loss  coefficient 
is  designated  £  .   Fig.  13  shows  the  associated  velocity  diagram  for 
a  turbine  stage. 

It  is  of  interest  here  to  summarize  the  assumptions  made  in  using 
the  one-dimensional  analysis,  as  they  will  have  a  large  bearing  on 
the  accuracy  of  the  results.   These  assumptions  are: 

1.  Adiabatic,  steady  flow,  so  that  the  absolute  and  relative  en- 
thalpies are  constant  respectively  along  any  streamline. 

2.  Cylindrical  stream  surfaces,  so  there  is  no  bending  of  the 
stream  surface,  or  radial  component  of  velocity. 

3.  Uniform  conditions  ahead  of  the  stator  so  there  is  no  total 
pressure  or  total  temperature  gradient  in  the  radial  direction. 

4.  Completely  radial  flow  into  the  plenum,  and  axial  flow  into 
the  stator  with  no  tangential  velocity  component. 

5.  Shear  stresses  are  negligible  on  the  annulus  and  shroud,  and 
symmetrical  wakes  occur  at  the  trailing  edges  of  the  blades. 
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The  validity  of  these  assumptions  will  be  discussed  later  and 
where  invalid  assumptions  were  found  to  exist  their  influence  on  the 
loss  coefficient  will  be  discussed. 

When  analysing  the  reduced  data,  nearly  all  variables  were  plotted 
as  functions  of  the  stator  discharge  pressure  ratio  P   /P  .   The 
following  discussion  will  show  why  P   /P   is  an  important  parameter. 
From  the  equation  of  momentum 

M  =  -^  RV 
gc    u 


and  using 


$  = 


WT  \Ag/gc  TT0 

P    A 
*T0   ths 


1  - 


n  -  1 


TO 


where 


n  = 


1  +  C   (Y  "  1) 
P 


and 


C„  " 


P   -  P 

TO    Tl 

P   -  P 
TO     1 


hence 


*  -  £ 


TO 


1  I  P, 


,  n,  y 


Thus 


M  = 


TO   ths 


8CN/V^ 


TO 


RV 


or 


PT0    \ef 


V^ 


TO 


TO 


1  \    P. 


n,  Y 
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Now 


or 


V   =  V   sin 
u 


V   =  C  T 
u       TO 


*  =  C(TT0  "  Tl) 


.  2 
sin  a 


n  -  1 


1  - 


.  2 
sin  a 


and 


sA 


=  =^ 


Hence 


TO 


M    =  f 


n  -  1 


1  - 


TO 


sin  a  =  f 


TO 


2  V  P, 


n,  YJ 


TO 


ref    3  V  P, 


n,  Y»  a, 


For  an  ideal  gas  y   is  constant,  and  for  small  pressure  ratios  n  is 
nearly  a  constant.   Thus,  if  it  is  assumed  that  the  stator  discharge 
angle  is  constant: 


M     c  -   f 


TO 


ref  4        P 


1. 


and 


and 


=   f 


TO 


11    P, 


=   f 


TO 


VTT0 


2    V    P, 


Using  the  loss  coefficient  as  defined  by  Vavra  [Ref.  8] 
2 


?s  =  1 


=  1  - 
1  -  1 


V   \  2    T-T      T/T   -T    /T 
1  \     1    lis    V    TO    lis'  TO 


V 


lis 


T   -  T 
TO    lis 


1  -  T    /T 

lis'  TO 


(ypTo)  n  -  jh'*^ 


Y  -  1 
Y 


1  - 


(VPTo) 


Y  ~  1 


=  f 


TO 


s  I  P 
1 


n,  Y 
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And,  consistent  with  the  above  assumptions  for  n  and  ,y, 

,  ■ ,  0 

B.   DESCRIPTION  OF  METHOD  1 

It  is  possible  to  solve  for  the  stator  discharge  properties  by 
satisfying  continuity  at  the  throat  plane  and  the  exit  plane,  provided 
certain  parameters  are  accurately  known  and  some  assumptions  are  made. 
At  the  exit  plane,  from  continuity: 

WT  =  pl  Val  Al 
From  the  perfect  gas  relationship 


P-i  = 


1   R  Tn 
8   1 


Thus 


WT  Rg  Tl 
V   =  — — s — -  (1) 

al    P  A  K   ' 

However  T  is  not  known  until  V   is  determined.   If  one  assumes  that 
the  same  pressure  and  temperature  acting  at  the  exit  also  exist  across 
the  throat  plane,  then  continuity  may  be  applied  at  the  throat: 

*T  =  pl  Vl  Aths 

Then  using  _ 

P  V 

Pl  =  R^T"  and   Tl  =  TT0  *  ~^~ 
8   1 

wT       w   R  T1     w„  R 
v   =  I =   T   8   L   =   T   g 

1   Pl  Aths     Pl  Aths    Pl  Aths 


or 


o    C  Pi  A^u  ^  R 

2   l_ths     m        T  g 

1     •   -        1    Pi  A^,     TO 
w^  R  1   ths 

T  8 
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Let 

C      -  1      ths 

2  WTEg 

Thus  „ 

Vl     +  C2  Vl   "  C   TI0 

And  completing  the  square: 

vi2  +  c2  vi  +   {-¥)     =  (vi  +    2)     ■  c  Tio  +  (4 


Hence 


$ 


vi  =V  ctto  +\~^l     ~     Hr-)  (2) 

and 

Vl 
T      =    T        -  ■ 

1  TO  C  (3) 

Once  V  1  is  determined  from  (1) , 
ai 

,/vA        •  /  A ,-.\ 

_1  I    il  l  _1   I    «-v.  <-.  \ 

and  V  ,  =  V,  sin  a, 
ul    1      1 

and  all  quantities  are  available  to  determine  C  . 

s 

It  is  important  to  examine  in  detail  the  assumptions  made  in 
arriving  at  c,     by  continuity  only,  which  are: 

1.  The  same  pressure  P  and  temperature  T  exist  at  the  throat 
and  exit  plane. 

2.  There  is  no  change  of  V   from  the  throat  to  the  exit  plane. 

3.  The  correct  pressure  distribution  for  P.  is  accurately  known. 

4.  The  areas  A  .   and  An  are  known  exactly. 

ths      1  J 

5.  There  is  one-dimensional  or  uniform  flow  at  both  the  throat 
and  exit  plane,  which  implies  no  boundary  layer  restriction. 

In  the  data  reduction,  P..  was  taken  as  the  mean  of  the  measured 

hub  and  tip  pressures ,  and  A  was  the  annulus  exit  area  reduced  by  a 

factor  K  to  account  for  the  trailing  edge  thickness. 
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C.   DESCRIPTION  OF  METHOD  2. 

In  this  method  both  continuity  and  momentum  are  applied,  where 
the  tangential  velocity  is  determined  from  the  momentum  equation, 
while  the  axial  velocity  is  determined  from  knowledge  of  the  axial 
forces  acting  on  the  stator  and  by  applying  continuity  at  the  stator 
exit.   From  angular  momentum: 
M  g 


V 


ul 


w  R 
T  ml 


(4) 


At  the  stator  exit  plane,  a  parabolic  pressure  distribution  is  assumed 
to  exist,  which  is  expressed  by  equation  12  of  Lenzini  [Ref.  6]  and  is 
repeated  here: 
P. 


hub 


lav 


(1   +   e.)    R   2      +  R    .      R,     .     -    (2   +   O    R   2L 
tip  tip      hub  hub 


(R 


tip 


w 


til 


(2   +  e)    R   2      -   R    .      R     ,     -    (1   +   e)    R,       2 
tip  tip     hub  Tiub 


(RtiP  "  *hub  > 


(5) 


The  value  of  e  determines  the  shape  of  the  parabola.   A  value  for  e 

is  first  assumed  and  with  this  first  approximation  for  P.   ,  the  force 

r  lav 

balance  equation  may  be  solved  for  V   , 


V 


alM 


w„ 


R  -  Pn    A 

lav  axs 


(6) 


where  R  ■  the  algebraic  sum  of  the  measured  forces  and  the  pressure 

forces  acting  on  the  stator.   Thus  V  ...  is  a  direct  function  of  the 
6  alM 

value  of  e,  or  the  shape  of  the  pressure  distribution.   Continuity  is 
then  written  for  the  exit  plane  as  in  Method  1,  employing  the  perfect 
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gas  relations,  to  obtain; 


ValC 
Let 


WTRg 
"P1A1 


wT  R 

C  =   T  « 

Ll   P,  A 


TO 


1   1 


WTRE 
P1A1 


"TO 


VaXC  +  Vul 


Then 


ValC  "  Cl  TT0 


aiC 


ul 


or 


2    C  2 

V  \      +   —  V    =  CTTn  -  V 
ale   C    ale     TO    ul 


and  completing  the  square: 


v  2  +  C_  v    + 
ale   C    alC 


m!  ■  c- 


+  2C. 


CT 


TO 


ul 


(m 


or 


V 


ale 


/CTI0  -  Vul  +  (4)2   "  (^) 


(7) 


Equations  (6)  and  (7)  are  compared  and  the  value  of  e  is  changed  appro- 
priately until  the  two  methods  produce  the  same  value  of  V  1 .   In  sum- 
mary then,  the  steps  are: 

1.   Assume  an  initial  value  of  e   and  determine  P,    from  equation 

lav 


(5). 


2.  Solve  for  V  ,  from  equation  (4). 

ul 

3.  Solve  for  V  ,..  from  equation  (6). 

alM       n 

4.  Solve  for  V  ,   from  equation  (7). 

5.  Compare  the  two  values  of  V   and  vary  e  until  the  two  values 


of  V  ..  become  equal, 
al 

6.   Compute  c,  . 
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The  assumptions  made  using  this  method  are  as  follows: 

1.  The  stator  torque  can  be  measured  accurately,  yielding  a  cor- 
rect value  of  V  .. 

ul 

2.  The  correct  pressure  distribution  for  P  is  known  accurately, 
or  conversely  the  iterated  value  of  e  will  yield  the  correct  pressure 
distribution. 

3.  The  correct  axial  force  is  being  recorded. 

4.  The  area  A  is  known  exactly. 

5.  There  is  one-dimensional  or  uniform  flow  in  the  nozzles,  so, 

no  boundary  layer  correction  factor  is  required. 

As  in  Method  1,  the  area  A.  was  reduced  by  the  same  factor  K  to 
•  1  t 

account  for  the  effect  of  trailing  edge  thickness. 

If  the  assumptions  for  both  methods  are  met,  then  either  method 

should  result  in  essentially  the  same  values  for  V  , ,  V  , ,  a.  and  c  , 

al   ul  s 

respectively. 
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.  V.   NATURE  OF  PROBLEM 

A  survey  of  previous  results  has  shown  that  £  has  ranged  from 
negative  values  to  as  high  as  0.30.   The  result  of  this  has  produced 
stator  or  rotor  efficiencies  of  over  one-hundred  percent,  but  the 
overall  machine  had  a  total-static  efficiency  in  the  order  of  eighty 
percent  at  the  design  point.   The  conclusion  reached  was  that  the 
determination  of  overall  performance  was  correct,  but  that  some  error 
existed,  either  in  the  theory  of  calculating  stator  loss  coefficients 
or  in  measuring  the  properties  necessary  to  arrive  at  the  loss  coef- 
ficient. 

Appendix  A  presents  several  methods  by  which  the  stator  discharge 
angle  and  stator  loss  coefficients  may  be  predicted.   In  addition,  a 
curve  fit  of  data  is  given  by  Vavra  [Ref.  8].   A  summary  of  these 
values  and  some  representative  ones  from  previous  work  by  Lenzini 
[Ref.  6]  are  presented  in  Table  I  below. 

TABLE  I 
Predicted  and  Experimental  Values  of  z, 


Predicted 
Values 

Source 

"s 

Ainley 
Traupel 
Soderberg 
Vavra 

0.0756 
0.1545 
0.1022 
0.1400 

Experimental 

Lenzini       '    -0.0100  to  0.3000 

1 
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Since  previous  work  by  Lenzini  [Ref.  6]  indicated  that  the  over- 
all total-static  turbine  efficiency  was  about  eighty  percent,  it  was 
concluded  that  many  of  the  above  values  were  too  high,  and  that  an 
average  value  for  the  loss  coefficient  c,     of  0.08  would  be  more  in 
order.   The  problem  then  became  one  of  determining  where  the  error 
in  computing  the  loss  coefficient  was  occuring  and  to  reduce  the  error 
as  much  as  possible. 

In  order  to  determine  the  error,  it  was  assumed  that  the  theory 

to  calculate  c,     was  correct  and  that  the  error  arose  from  invalid 
s 

assumptions  and  incorrect  measurements  of  some  variables.   If  a  value 
for  L,      is  assumed,  certain  parameters  can  be  calculated  and  plotted 
as  a  function  of  PTn/P '•,  >  and  then  actual  measurements  could  be  com- 
pared against  these  predicted  values.   Thus,  based  on  the  definition 

of  £  ,  several  important  relationships  could  be  developed: 

Y  -  1 
Tl  -  Tlis     (TT0  -  Vl2/C)  -  TT0  (V PT0)   " 

'S   "  It0   "  ""•  I       -  T       f?/*\     ^ 

TO  TO 


(vpto) 

v   -  1 
=  (VPT0)      ' 

__  (l   ~   Vl2/CTTo)   -   X     =        _   V] 


_-   1 
Let  X  = 


Hence 

17CTT0 

cs  :  ~T~T~  L  "(I  -  x) 


R 

v  g 

and   using  C      =  ' — -  -r2- 

P        Y   -    1   J 


h-  J~^t  0 - OC1  - x) 


With  the  chosen  value  of  r,     =   0.08  and  y  =  l.A, 

8 


28 


(8) 


v  1 

—    =   2.1450    (1  -  X)    2  (9) 

aT0 

and 

^-=1-^    (^)  (10) 

^0  \     TO  / 

From 
T       Pi      1      ths        PT0lPT0^Tl/        RgTI0        ^   ^  Ttq 


WIRgTT0  ;iA/8c   TI0        /PlY        lV"l 


P_„/\  T,/T_,A  a. 


PT0  Aths   >/«c  Rg  TT0  nA~  PT0  Aths  V  T0/  V    X     T0/  V  "" 


or 


TO/    \  V    TO    /  \aT0 


(ID 


WT  ^  P        H        V   - 

Then  M  =  —     R   ,    V  ths      T0     ml     ul 


g  ml      ul 


Sc    )/  Rg/gc   TT0 


or 


PT0  ref  ^  ths     ml      VaT0/  (12) 
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VI.   DESCRIPTION  OF  TESTS 

A  new  electrical  read-out  system  was  employed  for  the  strain- 
gauges  on  the  closure  plate.   It  was  decided  to  calibrate  both  the 
closure-plate  axial  force  and  torque  with  the  closure-plate  mounted 
on  the  stator,  rather  than  using  the  previous  method  of  calibrating 
away  from  the  TTTR  and  then  possibly  destroying  the  calibration  during 
the  installation  in  the  stator  casing. 

A  summary  of  the  operating  runs  is  listed  in  Table  II.  The  first 
series  of  runs  was  carried  out  at  pressure  ratios  of  1.5  and  2.0  over 
the  full  range  of  RPM  allowable  by  the  torque  capacity  of  the  air  dy- 
namometer. The  configuration  was  unhooded  with  rotor  axial  clearances 
of  0.250  in.  and  0.500  in.  Several  pressure  leaks  developed  on  nearly 
every  run,  and  only  Run  8  was  considered  sufficiently  trouble-free  to 
be  used  for  analysis. 

Once  the  high  loss  coefficients  were  confirmed  to  exist  from  the 
data  of  Run  8  and  reduction  by  Method  1  and  Method  2,  it  was  decided 
to  attempt  some  locked  rotor  tests  at  axial  clearances  of  0.250  in. 
and  0.500  in.  and  2.000  in.  to  determine  the  conditions  at  zero  RPM 
and  a  pressure  ratio  of  1.5.   Few  problems  were  incurred  during  this 
series  of  tests,  and  data  points  could  be  repeated  to  within  1  percent 
or  better.   During  the  last  run  however,  several  pressure  lines  were 
torn  loose,  which  required  extensive  repairs  before  any  further  runs 
could  be  made. 

During  the  runs  up  to  and  including  Run  14,  the  standard  method 
of  operation  had  been  to  leave  the  inlet  valve  of  Tank  1  full  open  and 
gradually  open  the  turbine  inlet  valve  to  increase  the  pressure  ratio. 
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TABLE  II 


Summary  of  Runs 


Run 

Configuration 

Purpose 

1-4 

ASME  Standard  Pipe  and  Orifice  Installed 

P    varied 
noz 

To  determine  nozzle 
calibration  coefficient. 

5 
6 

Rotor  removed 
Rotor  removed 

Closure  plate  axial  and 
torque  calibration 

7 

8 

9 

10 

11 

Operating,  PR  =  2.0,  AX  =  0.250" 
Operating,  PR  =  1.5,  AX  =  0.250" 
Operating,  PR  =  1.5,  AX  -  0.500" 
Operating,  PR  =  2.0,  AX  =  0.500" 
Operating,  PR  =  2.0,  AX  =  0.500" 

To  check  the  system 
after  modifications 
and  determine  operating 
data 

J 1 -■  ■  ' 

12  j Locked  Rotor,  AX  =  0.500" 

13  Locked  Rotor,  AX  =  2.000" 

14  (Locked  Rotor,  AX  =  0.250" 

To  examine  the  system 
operation  at  zero  RPM 
and  various  clearances. 

15   ^Operating,  PR  =  1.5,  AX  =  0.250" 

i 

! 

To  confirm  previous 
results 

16  JHoneycomb,  flexures  installed           |To  determine  the  ef- 

17  (Honeycomb,  flexures  removed             jfects  of  the  torque 

18  (Honeycomb,  flexures  installed           [flexures  and  various 

19  >Honeyoomb,  flexures  removed             j turbine  inlet  valve 

j openings . 

20  ! Nozzles  installed,  rotor  removed,       iTo  determine  the  ef- 
j                                     i 

flexures  removed               ifects  of  torque  flex- 

21  Nozzles  installed,  rotor  removed,        lures  and  inlet  valve 
i                                 i 

flexures  installed              ^openings  with  no  rotor 

i influence. 

! 
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As  a  result,  a  pressure  drop  in  the  order  of  50  in.  of  mercury  was 
occuring  across  the  turbine  inlet  valve.   Also,  it  was  noticed  that 
after  shutting  off  the  flow  to  the  stator,  there  were  residual  values 
of  stator  axial  force  and  torque  in  the  order  of  1.0  and  5.0  percent 
respectively,  which  gradually  dropped  to  zero  after  some  time. 

In  an  effort  to  examine  the  validity  of  the  assumption  of  no  in- 
let tangential  velocity  ahead  of  the  stator,  the  stator  was  removed 
and  the  honeycomb  flow  straightener  shown  in  Fig.-  9  was  installed. 
Also,  during  the  next  series  of  tests,  the  method  of  operation  was 
changed  and  a  further  modification  was  employed. '-For  Run  16,  the 
honeycomb  was  installed  and  the  turbine  valve  operated  as  before. 
Because  a  residual  torque  reading  was  again  observed,  it  was  decided 
to  remove  the  stator  torque  flexures  which  take  up  some  of  the  torque 
and  transmit  the  axial  force  to  the  cradle  assembly.   A  system  of 
roller  bearings  was  installed  between  the  cradle  assembly  and  stator 
plenum  casing  in  place  of  the  torque  flexures,  as  shown  in  Fig.  10. 
Also,  a  thrust  collar  was  fastened  on  the  stator  support  shaft  and 
allowed  to  press  against  the  support  bearing,  as  a  means  of  transmit- 
ting the  stator  axial  force  to  the  cradle  assembly.   The  next  series 
of  runs  was  performed  by  either  gradually  opening  the  turbine  valve 
to  obtain  the  desired  pressure  ratio,  or  fully  opening  the  valve  and 
adjusting  the  pressure  ratio  by  slowly  closing  a  dump  valve  in  the 
delivery  line  near  the  Allis-Chalmers  compressor.   It  was  possible 
to  obtain  similar  flow  rates  as  used  during  all  previous  tests,  but 
because  of  the  much  less  restricted  flow,  considerably  smaller  pres- 
sure ratios  were  required. 
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Prior  to  Run  16,  a  thermocouple  was  installed  on  the  torque 
capsule  support  assembly,  and  during  the  next  series  of  tests  the 
temperature  of  the  assembly  was  monitored.   Also,  the  discharge  tem- 
perature from  the  Allis-Chalmers  compressor  was  varied  from  90   to 
130   F  to  provide  a  sufficient  temperature  gradient  at  the  stator  in 
order  to  examine  the  effects  oi  differential  thermal  expansion.   Dur- 
ing Run  17,  the  flow  exit  angle  from  the  honeycomb  was  measured  using 
a  5-hole  probe,  to  find  out  whether  the  flow  leaving  the  honeycomb 
straightener  was  truly  axial  as  desired. 

Once  an  improved  method  of  operation  was  determined,  a  final 
series  of  tests  was  performed  with  stator  nozzles  reinstalled,  but 
the  rotor  removed  to  eliminate  any  influence  it  might  have.   The  same 
procedure  was  adopted  as  above,  with  the  flexures  installed  and  the 
turbine  valve  operated  either  fully  open  or  partly  open,  and  then  the 
flexures  were  removed  and  the  roller  bearing  supports  installed,  and 
the  valve  operated  either  fully  open  or  partly  open. 
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VII.   RESULTS  AND  ANALYSIS  OF  TESTS 

Fig.  14  shows  the  total-static  turbine  efficiency  as  a  function 

of  isentropic  head  coefficient.   It  is  seen  that  a  maximum  efficiency 

of  79.6  percent  occurs  at  a  k.   of  4.2.   Previous  results  for  these 

is 

operating  conditions  quoted  a  maximum  efficiency  of  83  percent  at  the 
same  k.  .   However,  the  efficiency  is  directly  proportional  to  the  dy- 
namometer torque  M  and  inversely  proportional  to  the  flow  rate  w 
for  given  inlet  conditions  and  speed  N.   Since  the  calculated  flow 
rate  is  now  considered  to  be  in  the  order  of  3.0  percent  higher  be- 
cause of  the  corrected  nozzle  coefficient,  this  would  account  for  the 
difference,  and  it  was  therefore  concluded  that  the  correct  results 
were  being  obtained. 

Fig.  15  shows  the  referred  stator  axial  force  as  a  function  or 
referred  RPM.   It  can  be  seen  that  there  is  some  non-linearity  when 
the  curve  is  continued  to  the  locked  rotor  condition,  which  can  be 
accounted  for  by  the  changing  pressure  distribution  around  the  shroud 
with  the  rotor  locked. 

Fig.  16  shows  the  referred  stator  torque  as  a  function  of  referred 
speed.   There  is  good  linearity  to  the  locked  rotor  condition,  indi- 
cating the  locked  rotor  torque  is  in  the  range  of  its  expected  value. 
When  the  referred  torque  is  plotted  as  a  function  of  stator  pressure 
ratio,  as  shown  in  Fig.  17,  there  is  a  difference  between  the  results 
of  Method  1  and  Method  2  and  a  considerable  difference  between  the 
actual  and  the  predicted  value  for  the  assumed  c,  =0.08.   The  dif- 
ference between  the  methods  is  due  to  the  resulting  higher  pressure 
of  P  by  Method  2.   It  can  be  concluded  however,  that  the  torque 
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measurements  at  a  pressure  ratio  of  1.45  are  about  8  percent  too 
high.   It  will  be  discussed  later  how  this  can  occur. 

Fig.  18  shows  the  flow  function  $  plotted  against  the  stator  pres- 
sure ratio.   The  curve  extends  well  to. the  locked  rotor  condition. 
The  reference  curve  shows  that  the  actual  $  is  more  in  error  at  the 
lower  than  at  the  higher  pressure  ratios  when  using  Method  1.   The 
difference  between  the  curves  of  Method  1  and  Method  2,  as  discussed 
previously,  is  due  to  the  higher  P  by  Method  2.   Although  $  is  a 

function  of  both  w^  and  A  .   for  fixed  inlet  conditions,  it  was  con- 

I      ths 

eluded  that  A  .   was  in  error.   Build-up  of  dirt  deposits  on  the 
ths 

blades  and  the  presence  of  an  increasingly  thick  boundary  layer  at 

higher  pressure  ratios  would  result  in  a  lower  effective  throat  area 

than  was  used. 

Fig.  19  shows  the  dimensionless  velocity  ratios  V  ,/a„n  and 

al   TO 

V  , /am_  as  functions  of  the  stator  pressure  ratio.   In  both  cases 
ul   TO 

the  curves  extend  smoothly  to  the  locked  rotor  condition,  and  there 

is  generally  good  agreement  between  the  two  data  reduction  methods. 

It  can  be  seen  that  both  V  -.  /a_„  and  V  -.  /a^n   axe.   lower  than  the  cal- 

ai   10      ul   TO 

culated  values  for  an  assumed  £   =0.08.   When  using  Method  1,  V  - 

is  determined  from  equation  (1) ,  and  a  low  value  of  V  .  could  be  due 

al 

to  a  high  value  of  pressure  P  or  the  reduction  factor  K  .   When  us- 
ing Method  2,  V   is  determined  by  varying  z   until  equations  (6)  and 
(7)  agree,  which  means  in  effect  that  P   is  varied.   Method  2  implic- 
itly assumes  that  V  ,  is  correct.   Examination  of  the  behaviour  of 

ul 

equations  (6)  and  (7)  when  reducing  the  data  revealed  that  equation 

(7)  generally  produced  a  lower  V   than  equation  (6),  and  its  value 

al 

changed  very  little  during  the  iteration.   The  strongest  variables 
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in  equation  (7)  are  P..  and  K  ,  and  a  low  value  of  V  ..  would  be  due 

to  high  values  of  P  and  K  ,  which  is  in  agreement  with  the  previous 

conclusion  for  Method  1.   The  magnitude  of  V  ,  has  little  effect  on 

ul 

V   in  equation  (7) .   Since  V  ,  by  equation  (6)  is  generally  higher 

than  by  equation  (7),  and  because  changing  P    causes  equation  (6) 

to  change  more  rapidly  than  equation  (7),  the  iteration  scheme  forces 

P-  higher  in  order  to  reduce  the  net  stator  force  and  thus  reduce  V  n 
1   &  al 

in  equation  (6)  .   The  increase  in  P..  causes  only  a  slight  reduction 

in  V   by  equation  (7) .   If  the  initial  value  for  Pn  with  e  =  0  were 
al  i 

too  high  then  the  above  undesirable  situation  would  occur,  and  the 

net  result  would  be  too  high  a  value  of  P..  and  too  low  a  value  of  V  .  . 

i  al 

A  desirable  condition  would  occur  if  the  initial  value  of  V  .  by 

al  J 

equation  (7)  were  higher  than  that  by  equation  (6) ,  since  the  iter- 
ation scheme  would  then  reduce  P..  causing  a  higher  resultant  V  1  . 
Previous  work  by  Commons  [Ref.  4]  indicated  that  the  stator  hub  pres- 
sure and  the  stator  tip  pressure  measured  with  the  arrangement  of 
Fig.  6  were  5  percent  higher  and  5  percent  lower,  respectively,  than 
the  pressured  obtained  by  flow  surveys  after  the  stator.   The  over- 
all effect  of  these  errors  would  produce  a  higher  initial  P.  than  the 
correct  value,  which  is  consistent  with  the  previous  conclusion  for 
the  low  value  of  V   .   Examination  of  the  results  from  Run  15  given 
in  Appendix  D  shows  that  the  hub  pressure  acting  on  the  hub  area  pro- 
duces the  largest  effect  on  the  net  force  in  equation  (6) .   As  a  con- 
sequence, a  slightly  low  value  of  P    would  result  in  the  undesirable 
condition  of  a  higher  value  of  V  1  by  equation  (6)  than  by  equation 
(7). 
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When  using  Method  1,  the  low  value  of  V  , /a„,A  in  relation  to  the 

°  ul   T(J 

calculated  curve  for  an  assumed  C   =0.08  would  be  due  to  the  low 

s 

values  of  V  n  and  Vn  caused  by  an  erroneously  high  pressure  P..  .   This 
al      1  1 

supports  the  previous  argument  where  it  was  stated  that  high  values 

of  P,  could  be  causing  the  errors.   Using  Method  2,  equation  (4)  shows 

that  the  low  value  of  V  ,  could  be  caused  by  a  low  stator  torque  read- 

ul 

ing.   However,  that  condition  would  conflict  with  the  conclusion  from 

Fig.  17  which  shows  the  torque  to  be  reading  about  8.0  percent  high. 

A  reduction  in  the  stator  torque  would  result  in  a  further  deviation 

of  V  ,  from  the  calculated  curve  for  C   =0.08,  and  an  increase  in 
ul  s 

the  mean  loss  coefficient.   This  situation  can  be  resolved  only  by 

concluding  that  the  angle  of  75.9   chosen  to  calculate  the  velocity 

components  from  equation  (9)  was  too  high.   A  reduction  in  the  chosen 

flow  angle  would  result  in  a  higher  V  n  and  lower  V  ,  for  a  given  V, . 

al  ul  1 

Fig.  20  shows  the  stator  loss  coefficients  calculated  by  the  two 
methods  as  a  function  of  stator  pressure  ratio.   It  can  be  seen  that 
the  losses  obtained  by  either  method  are  at  least  50  percent  higher 
than  the  predicted  value.   The  higher  loss  coefficient  by  Method  2 
is  due  to  the  higher  pressures  resulting  from  Method  2  and  the  fact 
that  constant  pressure  lines  diverge.   Fig.  20  also  shows  that  the 
loss  coefficients  determined  by  both  methods  decrease  slightly  as  the 
stator  pressure  ratio  and  the  flow  rate  increases.   Such  a  condition 
is  opposite  to  that  expected,  since  at  higher  pressure  ratios  an  in- 
creased possibility  of  separation  exists,  and  a  thicker  boundary  layer 
with  greater  secondary  losses  should  occur.   The  apparent  reduction  of 
the  loss  coefficient  as  the  stator  pressure  ratio  increases  is  likely 
due  to  an  unexpected  flow  condition  at  higher  flow  rates.   It  will  be 
shown  later  that  at  higher  flow  rates  a  changing  inlet  flow  condition 

does  occur. 

37 


Fig.  21  shows  the  referred  stator  axial  force  as  a  function  of 
stator  pressure  ratio  for  three  axial  clearances  of  0.250  in.  0.500 
in.  and  2.000  in.   It  can  be  seen  that  at  a  clearance  of  0.500  in.  a 
more  negative  axial  force  occurs  than  at  0.250  in.  or  2.000  in.   It 
is  suspected  that  the  presence  of  the  locked  rotor  at  various  axial 
clearances  causes  a  changing  pressure  distribution  on  the  shroud,  and 
thus  the  rotor  axial  clearance  has  a  large  effect  on  the  stator  axial 
force  and  the  shroud  pressure  distribution. 

Fig.  22  shows  the  referred  stator  torque  as  a  function  of  the 
stator  pressure  ratio.   The  same  reversal  of  trend  occurs  when  the 
axial  clearance  is  increased,  confirming  the  fact  that  the  presence 
of  the  rotor  influences  conditions  upstream.   The  locked  rotor  results 
shown  in  Fig.  22  at  the  three  different  axial  clearances  further  sup- 
ports the  previous  conclusion  from  Fig.  17,  namely  that  the  stator 
torque  is  reading  high. 

Fig.  23  shows  the  referred  flow  rate  as  a  function  of  stator  pres- 
sure ratio  for  the  three  axial  clearances  at  the  locked  rotor  config- 
uration. It  can  be  seen  that  there  is  good  agreement  at  each  of  the 
three  axial  clearances,  so  the  changing  torque  with  clearance  is  not 
due  to  a  different  flow  rate.  Fig.  23  also  shows  that  the  flow  rate 
approaches  a  maximum  as  the  pressure  ratio  nears  the  critical  value, 
but  that  choking  does  not  occur. 

Fig.  24  shows  the  flow  function  $  with  changing  stator  pressure 
ratio  for  the  locked  rotor  condition.   As  in  Fig.  23  there  is  no  ef- 
fect of  clearance  on  flow  rate.   However,  the  same  trend  of  a  lower 
flow  function  than  desirable  is  occuring,  which  confirms  the  results 

of  Fig.  18. 

38 


Fig.  25  shows  the  stator  discharge  pressure  ratio  as  a  function 
of  the  overall  pressure  ratio  P   /P  .   It  can  be  seen  there  is  a 
linear  relationship  between  stator  pressure  ratio  and  overall  pres- 
sure ratio.   Fig.  25  also  shows  the  increasing  difference  between 
Method  1  and  Method  2  as  the  overall  pressure  ratio  is  increased. 

Fig.  26  shows  the  referred  velocities  as  a  function  of  stator 
pressure  ratio  for  the  locked  rotor  configuration.   The  same  trend 
in  V  1  occurs  with  the  locked  rotor  results  as  was  noticed  in  Fig.  19 
when  the  rotor  was  unlocked,  confirming  the  conclusion  that  some  er- 
ror exists  in  the  pressure  P  and  the  reduction  factor  K  .   However, 
Method  2  results  in  a  tangential  velocity  component  which  is  greater 

than  the  calculated  value  for  the  assumed  £   =0.08.   It  can  there- 

s 

fore  be  concluded  that  the  increased  value  of  V  ..  /am_  is  directly 

ul   TO 

due  to  an  error  in  the  stator  torque  measurement  and  that  the  high 
torque  reading  is  caused  by  some  combination  of  temperature  effect 
on  the  assembly  and  an  inlet  swirl  at  the  stator  plenum. 

Fig.  27  shows  the  stator  loss  coefficient  as  a  function  of  stator 
pressure  ratio  for  the  locked  rotor  condition.   Equation  (8)  gives 
the  relationship  between  the  pressure  ratio,  exit  velocity  and  stator 
loss  coefficient,  and  it  can  be  seen  that  for  a  given  pressure  ratio 
a  low  velocity  will  result  in  a  high  loss  coefficient,  and  conversely, 
if  the  velocity  were  too  high,  the  loss  coefficient  could  become  neg- 
ative.  When  using  Method  1  the  loss  coefficient  is  a  unique  function 
of  the  throat  pressure  and  area  where  the  throat  pressure  was  assumed 
to  be  the  same  as  that  existing  at  the  exit  plane.   Behaviour  of  equa- 
tion (2)  shows  that  a  low  velocity  will  result  if  a  combination  of 
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the  pressure  and  throat  area  is  too  large.   It  has  been  previously 
shown  that  the  throat  area  appears  to  be  too  large  and  Figures  20  and 
27  show  that  at  low  pressure  ratios  the  throat  pressure  must  be  high 
also.   Appendix  C  gives  a  predicted  pressure  distribution  at  the  exit 
plane,  and  it  can  be  seen  that  a  parabolic,  rather  than  a  linear  var- 
iation would  be  more  in  order,  resulting  in  a  slightly  higher  P..  than 
the  linear  variation  from  hub  to  tip  would  predict.   From  this  it  may 
be  concluded  that  assumption  1  of  Section  4.B,  namely  that  equal  pres- 
sures exist  at  the  throat  and  exit  plane  would  be  in  error.   As  the 
stator  discharge  pressure  ratio  increases  the  loss  coefficient  de- 
creases somewhat,  indicating  that  the  pressure  at  the  throat  and  exit 
plane  are  more  nearly  equal.   When  using  Method  2,  an  erroneously  high 
stator  torque  reading  will  produce  a  high  tangential  velocity  compo- 
nent resulting  in  a  higher  net  velocity  than  should  exist  for  a  given 
pressure  ratio.   It  should  be  noticed  that  the  loss  coefficient  by 
Method  2  increases  with  increasing  pressure  ratio  even  though  the 
stator  torque  is  high.   This  indicates  that  the  error  in  hub  pressure 

P,  ,  and  exit  pressure  P.  have  a  stronger  effect  than  does  the  stator 
hub  1 

moment  and  it  was  concluded  earlier  from  the  analysis  of  Fig.  19  that 
P  was  indeed  too  high. 

Fig.  28  shows  the  measured  temperature  of  the  torque  capsule  sup- 
port assembly  as  the  flow  rate  is  increased,  and  Fig.  29  shows  the 
results  of  the  temperature  analysis  on  the  cradle  assembly.   It  can 
be  seen  that  as  the  cradle  assembly  heats  up,  a  positive  error  occurs 
in  the  stator  torque  reading,  indicating  that  the  reading  becomes  too 
high  as  the  inlet  temperature  increases.   The  reason  for  the  increasing 
torque  error  is  the  unequal  expansion  of  the  capsule  support  assembly 
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and  the  connectors  between  the  point  of  action  of  the  torque  arm  and 
the  upper  end  of  capsule  support.   The  stator  torque  flexures  pre- 
vent free  rotation  of  the  stator  as  the  capsule  assembly  grows  due  to 
thermal  expansion.   The  two  steel  rods  and  the  Invar  ring  of  the  torque 
capsule  do  not  expand  by  an  equal  amount  in  the  opposite  direction  to 
cancel  the  thermal  expansion  of  the  capsule  support.   The  differential 
thermal  expansion  causes  an  additional  tensile  force  on  the  capsule 
and  a  higher  apparent  torque  reading.   Additional  analysis  showed  that 
it  would  be  impractical  to  compensate  for  this  thermal  expansion  by 
making  the  support  assembly  of  the  same  type  of  material  as  the  capsule 
ring  and  steel  attachments,  since  the  compensation  would  be  possible 
only  for  a  fixed  temperature.   In  addition  it  would  be  impractical  to 
thermally  shield  or  cool  the  assembly  as  was  done  with  the  dynamometer 
capsule,  since  this  set-up  would  present  severe  complications  if  the 
turbine  were  operated  inside  the  closed  hood.   Thus  the  only  practi- 
cal way  to  eliminate  the  thermal  effect  on  the  capsule  support  assembly 
was  to  give  the  stator  freedom  to  rotate  so  that  no  additional  tension 
could  be  applied  to  the  capsule.   A  solution  was  accomplished  by  re- 
moving the  torque  flexures  and  installing  roller  bearing  supports  to 
maintain  the  stator  alignment,  and  a  thrust  collar  was  fastened  to  the 
stator  support  shaft  and  allowed  to  press  against  the  support  bearing 
as  a  means  of  transmitting  the  stator  axial  force  to  the  cradle  assem- 
bly. 

Fig.  30  shows  the  referred  stator  torques  obtained  with  the  honey- 
comb flow  straightener  replacing  the  guide  vanes.   Test  data  are  shown 
at  different  openings  of  the  turbine  inlet  valve  and  the  tests  were 
carried  out  with  the  torque  flexures  installed  or  removed.   If  the 
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flow  after  the  honeycomb  were  axial,  and  if  the  flow  entered  the  sta- 
tor  assembly  with  no  whirl  components,  the  stator  torque  would  be  zero 
for  any  pressure  ratio.   Fig.  30  shows  that  the  stator  torque  was  only 
zero  at  zero  flow  rate,  and  increased  with  increasing  pressure  ratio. 
From  the  results  of  Fig.  30  it  can  be  concluded  that  temperature  ef- 
fects caused  by  the  presence  of  the  torque  flexures  and  throttling  the 
flow  by  partly  closing  the  turbine  inlet  valve  influence  the  stator 
torque  readings.   With  the  turbine  valve  full  open  less  torque  occurs 
when  the  flexures  are  removed,  indicating  that  the  temperature  effect 
due  to  the  flexures  has  been  reduced.   When  the  pressure  ratio  is  set 
with  the  turbine  valve  partly  closed  a  considerably  larger  torque 
occurs  and  it  is  again  higher  if  the  torque  flexures  are  installed. 
Since  an  increasing  torque  with  increasing  pressure  ratio  still  occurs 
when  the  turbine  valve  is  full  open  and  the  torque  flexures  are  re- 
moved, it  can  be  concluded  that  either  the  honeycomb  is  not  truly 
axial  or  the  flow  still  enters  the  stator  plenum  with  a  small  whirl 
component . 

Fig.  31  shows  the  flow  angles  measured  after  the  honeycomb 
straightener .   It  can  be  seen  that  the  velocity  after  the  honeycomb 
has  a  tangential  component  which  increases  with  increasing  flow  rate. 
From  Figures  30  and  31  it  can  be  concluded  that  the  partly  opened  tur- 
bine valve  causes  whirl  components  in  the  delivery  pipe,  and  the  flow 
enters  the  stator  assembly  with  a  tangential  velocity.   Hence,  assump- 
tion 4  of  Section  4. a,  namely,  that  the  flow  into  the  plenum  is  radial 
and  no  whirl  component  exists  ahead  of  the  stator  is  not  satisfied 
with  the  present  installation.   Because  the  honeycomb  has  openings  of 
about  1/8  inch  it  is  apparently  not  capable  of  completely  turning  the 
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flow  into  the  axial  direction.   It  was  found  also  that  the  stator 
blockage  factor  E,   through  the  honeycomb  is  about  0.72,  where  g  rep- 
resents the  ratio  of  the  actual  flow  area  and  that  necessary  to  pass 
the  flow  if  the  expansion  were  isentropic.   Because  of  the  low  valves 
of  E,  it  can  be  concluded  the/e  must  exist  large  whirl  components  ahead 
of  the  honeycomb  causing  flow  separation  at  the  inlet  and  large  losses, 
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VIII.   CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  preceding  analysis,  it  can  be  concluded  that  the  stator 
torque  reading  is  about  8  percent  too  high,  and  that  this  error  is 
due  to  a  combination  of  the  temperature  effect  on  the  capsule  support 
assembly  and  an  inlet  swirl  caused  by  the  turbine  inlet  valve  and 
piping.   It  is  recommended  that  the  turbine  be  operated  with  the  tur- 
bine inlet  valve  full  open  at  all  times,  and  that  the  pressure  ratio 
be  set  by  closing  a  dump  valve  near  the  air  supply. 

During  the  last  series  of  runs  with  the  stator  nozzles  installed 
and  the  rotor  removed,  inconsistent  and  residual  values  of  stator 
axial  force  occured  when  the  torque  flexures  were  removed  and  the 
roller  bearings  installed.   This  was  due  to  the  mechanical  friction 
inherent  in  the  present  system.   It  is  felt  that  a  better  system 
should  be  designed  which  would  allow  the  stator  to  rotate  and  move 
axially  without  restraints.   This  would  make  the  stator  a  truly  free- 
floating  assembly  and  eliminate  the  influence  of  temperature  on  the 
measurement  of  the  stator  torque.   Such  a  system  would  require  the 
use  of  ball-bearings  mounted  in  such  a  way  that  they  would  allow  both 
rotation  and  axial  movement  of  the  stator  and  thermal  expansion  or 
contraction  of  the  casing,  but  maintain  the  stator  axial  alignment. 

The  analysis  of  the  cradle  assembly  showed  also  that  thermal  ef- 
fects occur  in  the  horizontal  direction  which  affect  the  capsule  meas- 
uring the  stator  axial  force.   This  problem  could  be  eliminated  by 
constructing  a  cradle  support  assembly  similar  to  that  used  for  the 
torque  measurement  but  made  of  aluminum.   Proper  installation  of  this 
cradle  would  allow  its  growth  or  contraction  to  be  in  the  same 
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direction  as  the  horizontal  members  of  the  stator  cradle,  thus  com- 
pensating for  thermal  effects  on  these  members. 

It  is  suggested  that  a  honeycomb  flow  straightener  and  series  of 
screens  be  installed  in  the  delivery  pipe  just  ahead  of  the  inlet  to 
the  stator  plenum  chamber,  ana  that  the  conical  screen  presently  in- 
stalled in  the  plenum  chamber  be  removed.   This  should  produce  a  more 
uniform  flow  into  the  plenum  and  eliminate  the  tangential  component, 
while  removal  of  the  conical  screen  should  eliminate  some  of  the 
existing  total  pressure  gradient  at  the  entrance  to  the  stator  noz- 
zles . 

From  the  analysis  of  the  pressure  effects  on  the  loss  coefficient 
it  can  be  concluded  that  extensive  flow  surveys  should  be  done  at  dif- 
ferent operating  conditions  to  accurately  determine  the  pressure  dis- 
tribution at  the  stator  discharge.   Means  must  then  be  found  to  ensure 
that  accurate  hub  and  tip  pressures  can  be  recorded  since  they  have  a 
large  effect  on  the  loss  coefficient.   It  was  shown  that  the  assump- 
tion of  stator  exit  pressure  also  acting  at  the  throat  plane  was  er- 
roneous and  the  assumption  could  be  eliminated  if  the  throat  hub  and 
tip  pressure  taps  were  accurately  located  at  the  throat  plane.   At 
present  there  is  up  to  a  5  percent  error  between  the  four  sets  of 
taps.   The  flow  surveys  must  be  carried  out  with  probes  which  prop- 
erly compensate  for  Mach  number  effect,  depth  of  immersion  and  pitch 
angle.   Very  useful  information  could  be  gained  by  using  a  hot  wire 
anemometer  probe  to  accurately  determine  velocity  profiles  and  mag- 
nitudes in  the  stator  area  at  different  operating  conditions. 

Results  of  the  prediction  of  closure  plate  torque  and  shroud  torque 
given  in  Appendix  B  show  these  values  never  exceed  0.1  percent  of  the 
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stator  torque  reading.   This  amount  is  considered  negligible  and 
thus  the  closure  plate  torque  webs,  strain  gauges  and  wires  could 
be  removed.   This  would  eliminate  the  congestion  in  the  hub  area 
and  allow  the  closure  plate  greater  freedom  of  movement,  since  some 
friction  might  be  caused  by  the  presently  installed  equipment.   The 
predicted  torque  as  a  function  of  pressure  ratio  and  speed  agreed 
very  well  with  the  strain  gauge  readings  and  so  it 'is  concluded  the 
closure  plate  and  shroud  torques  could  be  accounted  for  during  the 
data  reduction  by  using  an  empirical  formula.   Results  of  the  forces 
on  the  closure  plate  due  to  the  pressure  distribution  on  the  face 
agreed  to  within  1.0  percent  of  the  strain  gauge  reading.   Consid- 
erable effort  was  required  to  numerically  integrate  the  pressure 
over  the  face  area,  and  since  measurement  of  the  closure  plate  axial 
force  by  the  strain  gauges  was  much  quicker  and  produced  comparable 
results,  it  is  recommended  the  closure  plate  pressure  lines  be  re- 
moved. 

It  can  be  concluded  that  several  of  the  assumptions  utilized  to 
solve  the  equations  for  a  one-dimensional  analysis  were  not  valid, 
and  that  ignoring  their  effects  can  produce  very  poor  results.   Appen- 
dix C  shows  that  cylindrical  stream  surfaces  cannot  exist,  and  a 
radial  component  of  velocity  must  be  taken  into  account.   Flow  surveys 
made  by  Commons  [Ref.  4]  show  that  there  is  a  variation  of  the  total 
pressure  at  the  stator  entrance  and  thus  a  radial  entropy  gradient 
exists.   The  results  of  this  investigation  have  shown  that  the  flow 
is  not  completely  radial  into  the  plenum  or  axial  into  the  stator 
and  that  a  tangential  component  of  velocity  does  exist.   However, 
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these  unsatisfactory  conditions  can  be  reduced  by  incorporating  the 
equipment  modifications  suggested. 
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FIGURE  6 

TURBINE   AND  SHROUD     DETAILS 
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FIGURE  12 

THERMODYNAMIC     PROCESS      FOR  AN  AXIAL 
TURBOMACHINE 
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APPENDIX  A 

PREDICTION  OF  STATOR  DISCHARGE  ANGLE 
AND  STATOR  LOSS  COEFFICIENT 

1.  Stator  Discharge  Angle 

The  prediction  of  the  stator  discharge  angle  uses  methods  described 
by  Vavra  [Ref.  8]  and  Ainley  [Ref.  2].   In  these  methods,  it  is  assumed 
that  the  discharge  angle  is  influenced  only  by  blade  geometry  and  Mach 
number,  and  that  flow  incidence  angles  have  no  effect. 

For  values  of  M  _  0.5,  a  formula  given  by  Vavra  [Ref.  8]  is  used 
for  a  first  approximation,  and  the  effects  of  blade  geometry  and  Mach 
number  are  then  accounted  for  by  using  the  method  described  by  Ainley 
[Ref.  2].  A  linear  variation  of  a  is  then  assumed  for  values  of  M- 
between  M  =0.5  and  M  =  1.0.  The  results  of  these  calculations  for 
the  blade  geometry  shown  in  Fig.  7  are  listed  in  Table  III. 

2.  Stator  Loss  Coefficient  By  Method  of  Ainley 

The  prediction  of  stator  loss  coefficient  uses  the  method  described 
by  Ainley  [Ref.  2].   In  this  method  it  is  assumed  that  the  loss  coef- 
ficient is  independent  of  Mach  number  and  is  influenced  only  by  blade 
geometry  and  flow  incidence.   The  losses  are  divided  into  profile 
losses  Y  which  include  mixing  losses,  secondary  losses  Y  ,  and  tip 


clearance  losses  Y,  ,  where 
k 


P(i  =  0)  ^ 
(t/c\  al/a2 
and 


Y  + 

p(al  =  0) 


2  r~ 
Y 


-  y 
p  (al  =  a2)     p (al  =  0) 
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Y  +  Y  = 
s    K 


A  +  B 


(5) 


C  \  2   /  cos   ( 

s/cy 

cos   a 


m 


Values  of  Y  =  f(s/C)  are  given  in  Fig.  4  of  Ref.  2  and  values 


for 


X   =  f 


(VO 


_/i  +  i.d./o.dA 

are  given  in  Fig.  8.   A  relation  given  by  Vavra  [Ref.  9]  is  then  used 
to. relate  the  stator  loss  coefficient  to  the  pressure  loss  coefficient, 
The  results  of  these  calculations  are  listed  in  Table  IV. 
3.   Loss  Coefficient  By  Method  of  Traupel 

The  method  of  Traupel   is  one  of  the  most  complete  methods  used 
and  takes  into  account  a  large  number  of  variables  considered  to  have 
an  influence  on  the  loss  coefficient. 

p    w    r    zus 

where    z,  =   Total  loss  coefficient 

L  =  Profile  loss  coefficient 

P 

5    =  Wall  loss  coefficient 
w 

C    =  Secondary  loss  coefficient 

£    =  Damping  wire  loss  coefficient 
zus 

The  profile  loss  coefficient  is  given  by: 


'pO  '  Xp0 


•  XM  *  Xs  +  Sn  +  CF 


where    £  rt  =  Basic  loss  coefficient 
pO 

•  X  _  =  Reynolds  number  correction 


Traupel,  W.,  "Thermische  Turbomaschinen ,  Erster  Band,"  Springer 
Verlag,  Berlin/Goettingen,  Heidelber,  1958,  pg.  269-298. 
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y„  =  Mach  number  correction 
M 

X  =  Trailing  edge  thicftness  correction 

z,     =  Mixing  loss  correction 
m 

£_  =  Correction  for  straight  blades 

r 

The  wall  loss  coefficient,  which  is  due  to  the  friction  loss  at  the 

annulus  of  the  hub  and  tip,  is  given  by: 

s.  sin  a 
Cw  =  So  *  Xp0     h 

The  secondary  loss  coefficient  is  given  by: 

^r  =  *£  ^rO  +  Cs 

where     <K   =  Aspect  ratio  correction 

L    _  =  Basic  secondary  loss 
rO 

z,        =  Tip  leakage  loss 

A  summary  of  the  calculations  is  listed  in  Table  V. 
4.   Loss  Coefficient  Method  By  Soderberg 

The  method  of  Soderberg  is  taken  from  a  paper  by  Horlock. 
Soderberg  bases  the  loss  on  blade  geometry,  aspect  ratio,  thickness 
ratio,  and  Reynolds  Number,  and  then  relates  the  loss  coefficient  to 
the  gas  deflection  angle  through  the  aspect  ratio  only.   He  ignores 
such  parameters  as  Mach  number  and  certain  non-dimensional  gas  prop- 
erties. 

The  results  of  a  calculation  for  the  converging  TTTR  nozzles  is 
listed  in  Table  VI. 


Horlock,  J.  H.,  "Losses  and  Efficiencies  in  Axial  Flow  Turbines," 
Ins.  J.  Mech.  Sci. ,  Vol.  2,  1960,  pg.  59-60. 
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TABLE  III 

Prediction  of  Stator  Discharge  Angle 

s  =  0.8594  in 

a  =  0.2050  in 

t  =  0.0240  in 

a/s  =  0.2385  in 

t/s  =  0.0279  in 

For  M  £  0.5: 

*      "I  (   a/A  [Ref.  8] 

a  =  cos    r"K~) 

where  K  =  1  -  ^r  [  -  100  J   - 

t     io3  Vs   /   s 

=  0.9809 
Hence  a*  =  cos"1  (0.24318)  =  75.92° 


* 

a  =  a  + 

<\) 

where  e  = 

? 

8Z 

=  4.5338 

4(s/e)  = 

0.76° 

Hence  a  = 

=  76.68° 

For  M  = 

1.0 

[Ref.  2] 


a  (M  =  1.0)  -  cos"1  [ ~\=   cos  1    (0.27152)  =  74.25° 

s  -  t/cos  a  J 

For  M   (average)  =0.65  a   (average)  =  75.90 


82 


al : 

«  (T 

a2  = 

-75.90° 

t  = 

0 

c  = 

1 

4313 

s  ■ 

0 

8594 

s/c 

= 

0.6004 

TABLE  IV 


Loss  Coefficient  by  Method  of  Ainley 


Y  ,      .  -  0.0490 
p (al  =  0) 


Y  +  Y„  = 
s    K 


A  +  B  I  = 


i) 


s/C 


2 

cos   a, 

cos   a 
m 


Fig.  4 


tan  a  +  tan  a_ 

where  tan  a  =  = =  -1.9906 

m  2 


a  =  -63.33 

m 

B  =  0.5   For  non-shrouded  wheels 

°L 

— 27r  =   2    (tan  a,    -    tan   a,,)    cos    a      =    3.5720 

s/C  1  2m 


A   =   f 


(vaiV 

(1   +  I.D./O.D.) 


A2        AQ2   cos    a2 


Ai     Aoi  cos  ai 


ttD   _  h„    cos   a_ 

m2      2 l_ 

rrD   .,  h.,    cos   a- 

ml  1               1 


hl   =  V      °ml   =   °m2 


Hence   A   =    f 


2         ,        2 
cos      a9/cos      a, 

C1  +  "hubi^tipi). 


=   f    (0.0321)    =  0.0056 


Fig.    8 
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K  =  0 

2 
Y     +  Y     =  0.0056    (3.5720)2  cos      (-75.90)      =  00471 

S  cosJ    (-63.33) 


Hency  Y  =  Y     +Y     +  Y,    =  0.0961 
p  s  k 

t/s   =  0.0279 

Trailing  edge  correction  factor  F  =  1.02 

Hence  Y  =  1.02  (0.0961)  =  0.0980 


,fc 


Y  -  1 

Y 


1  +  Y 


Y(VPT0)_ 


-  1 


Y  -  1 


.VPT0 


-  1 


[Ref.  9] 

page  99 


With  M1  (average)  =  0.65, 


TO 


=  0.7055  (Run  15) 
average 


C   =  0.0756 
s 
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TABLE  V 


Loss  Coefficient  By  Method  of  Traupel 

C  „  =  f(a.   a.)  =  0.0040  Fig.  B3 

pO       U ,   1 

^po  =  fi(Re>  c"> 

Re  .  V  .  ?50  x  1.4313/12  =  5>72  x  1Q5 
V       1.564  x  10 

K  =  0 

Hence  x  n   =   0-92  Fig.  B5 

pO 

ik.  =  f(Mj  =  0.075  for  M.  =  0.65  Fig.  B4 

Ml  1 

^  =  f(l-K) 

1  -  K  =  1  -  (s — -\  =  1  -   /0.8594  -  Q,024° )  =  0.2392 

\    C°S  al)  \        cos  75.9°/ 


Hence  ijj  =1.60  Fig.  B6 

i,     =  f (1  -  K)  =  0.0700  Fig.  B6 

Cp  =  f(h/Dm) 

h/D  .  !Lm  =o.0814 
m   8.48 


Hence  £_  =  0.0250  Fig.  B7 

r 

Thus  c  -  0.0040  x  0.92  x  0.075  x  1.60  +  0.0700  +  0.0250 

P 

or  5  -  0.0954 

P 

r  n    nn/n         n    ao         0.8594    x   sin    75.9    _    „    _. . . 

C      =   0.0040   x  0.92   x  _    ,-.. =   0.0044 

w  0.690 

V*Po =  f(c/h) 
c/h  =  ¥li  = 2-078 
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Hence   ty  Jty        =   1.70  ELg.    B9 

and  ty     =   1.70    (0.92)   =  1.563 


?r0=  f 


*) 


v 

1     7  9^ 

-^  =  j^g  =  1.306  From  point  8  Run  15 


Hence  X,   _  =  0.0350  Fig.  B8 

rO  ° 


5   =  f(h  /h)  =  0  Not  applicable  to  stator 


Hence  z,     =  1.563  x  0.0350  =  0.0547 
r 

£    =0  Not  applicable  to  stator 

zus 


Hence  5  =  0.0954  +  0.0044  +  0.0547  =  0.1545 
s 
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TABLE  VI 


Loss  Coefficient  By  Method  of  Soderberg 


c1 

=  (  1  +  K  ) 

ao 

=  0 

al 

=  75.90° 

b  = 

0.975 

h  = 

0.690 

t/c 

=  0.248 

0.975 


+  0,  ::<{]-! 


Blade  axial  chord 
Blade  height 
Thickness /chord 


i      =  f(Aa)  =  0.0700 
Hence  £'  =  0.1580 


Fig.    5 


And   z,     = 
s 


1  +  5 


Y    -    1 
Y 


L1  +  I  CP1/PT0)J 


-  1 


VPT0 


Y    -    1 

Y  "    1 


With  M      (average)    =   0.65, 


=  0.7055    (Run  15) 


TO/       average 


%     =  0.1022 
s 
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APPENDIX  B 

PREDICTION  OF  TORQUE  ON  CLOSURE 
PLATE  AND  SHROUD 
1.   Torque  on  Closure  Plate 

Because  of  the  proximity  of  the  rotor  and  its  high  speed  rotation, 
a  core  of  rotating  fluid  will  be  set  up  between  the  rotor  and  closure 
plate,  and  friction  will  thus  transmit  some  torque  to  the  closure 
plate. 

Schlichting  '  gives  a  formula  for  the  torque  transmitted  to  a 

stationary  flat  plate  by  a  rotating  parallel  plate,  assuming  the  core 

of  fluid  rotates  at  one-half  the  speed  of  the  wheel. 

For  turbulent  flow: 

M  =  0.037  pa)  R   [  — -  \ 
V  ojR  I 


Using  values  of  p  and  v  for  air  at  520  R  and  one  atmosphere, 


with  R  =  3.87  in 

M 


-0.037x0.002378(11-)   (W 


173  x  10  6  x  60  x  144 
ttN  x  3.87 


1/5 


This  results  in  a  moment  of  0.007  ft-lbs  at  10,000  RPM  and  0.023  ft-lbs 

at  20,000  RPM. 

2 .   Torque  on  Shroud 

A  torque  will  be  transmitted  to  the  rotor  shroud  by  the  tangen- 
tial component  of  V1  and  by  the  scrubbing  of  rotor  at  the  tip.   A 


Schlichting,  H.,  "Boundary  Layer  Theory,"  pg.  607 
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g 

formula  given  by  Schlichting   for  the  torque  transmitted  between  two 
rotating  cylinders  is: 


u(r)  = 


du 


r   —  r 
2     1 


dr     2 


r2  "  rl   «-■ 


2   2 

2       2\    rl   r2 

a)2r2   -  to1r1  ) - ^  u2   -  u 


2   2 
r2 


(  W2  "  Wl)_ 


(w2r22  '  Yl2)   +   1r2   2  (U2  "  Wl)_ 


T  =  U 


du 
dr 


With  oj  =  0  and  r  =  r  , 


M  =  T  A2r2  =  m 


©C->0 


or  M  =  4iryh 


2   2 

rl   r2  "l 

2     2 
r   —  r 
2     1 


As  before,  with  the  core  assumed  to  rotate  at  one-half  the  rotor  speed, 
and  the  values  for  air  at  standard  conditions,  the  moment  transmitted 
is  0.018  ft-lbs  at  10,000  RPM  and  0.037  ft-lbs  at  20,000  RPM. 


Schlichting,  H.,  "Boundary  Layer  Theory,"  pg.  81 
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APPENDIX  C 


PREDICTION  OF  STATOR  DISCHARGE 
PRESSURE  DISTRIBUTION 
The  equation  of  motion  for  absolute  flow  through  a  stator,  as 
given  by  Vavra  [Ref.  9]  is: 
3V 


+  VH  =  V  x  (V  x  V)  +  T  Vs  +  f  - 
3t  f 


(CI) 


If  steady,  frictionless  flow  is  assumed,  the  equation  reduces  to: 

VH  =  V  x  (V  x  V)  +  T  Vs  (C.2) 

Each  variable  in  the  above  equation  may  then  be  expressed  in  an  axisy- 
mmetric  coordinate  system,  where 


Xm  3m    R   36    1n   3n 


.   —  da    e  as  .  -   3s  j 

m  3m   R   39    n   3n  /  . 


and  V  x  (V  x  V)  =  i  A  — 

m    K 


3(RVj    3V 
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36     3n 
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Equating  components: 

-  .   3H  _  f  !i  P(RV     ^m 

Xm*   3m  "  ;  R     3m         36 


-  V 


3V 
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+  V   K 
3n    mm 
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90 


3V 
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+  V   K 


n  [  8m    n  n_| 
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_9s 
9m 


(C.3) 


v 

-  ,  1  dU   m   __m 
e:  R  36  "  R 


3V    9(RV.) 
m 


96 


9m 


3V    9(RV  ) 
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90 


9n 


1  ds_ 

~  R   90 


(C.4) 


-  .   _9H 
U    9n 


-  V 
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9(RVQ)    9V 


3n 


90 


9V 


r-a  +  V  K 
dm    n  n 


+  V 


m 


9V 


m 


+  V  K 
9n    mm 


+ 


_9s 
9n 


(C.5) 


where     m  =  streamline  coordinate 

n  =  coordinate  normal  to  streamline 


K  =  curvature  of  streamline 
m 

K  =  curvature  of  normal 
n 


If  axisymmetric  flow  is  assumed,  all  terms  with  9/90  vanish  and  since 

the  equation  is  written  along  a  streamline  where  V  =0 

n 


i    _ffi  . .  0  -  Za    3<RV 

R         90 


9m 
9(RVJ 


Since  V   is  not  zero,  — - — 

m  dm 


=  0 


(C.6) 


(C7) 


3H    3«?  3H 

Then  —  =  -r—  and  if  adiabatic  flow  is  assumed  then  -r—  =  0  means  H  is 

dm    dm  dm 

constant  along  a  streamline.   In  addition,  —  =  0  means  the  flow  is 

dm 

isentropic  along  a  streamline.   A  consequence  of  these  last  two  state- 
ments results  in  isentropic  flow  or  V  x  (V  x  V)  =0.   The  remaining 
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equation  is  then: 


+  V„  -*=  +  V„  K  +tt  (C.8) 


3n   R    3n        m   3n     m   m   3n 


or 

a  /  V  2\     9      Vfi  9<RV  sh       ^o 

.J.   JH-U  V  2  K  +  /  0-  -  M  +  T  M  .  0 

3n  \   2  /   m   m   R  n  3n       3n 


V  +  vm2. 

Employing  the  relation  H  =  C  T  +  \ J  the  following  important 


( 


equation  results: 


dn    \     m     J  mVm   C    3ny 


Vfl   3(RV  )  3T 

+  2  — 2  C   — - 

•  R    3n       L     p   3n 


P 
+  (2TT"V  )  ^=°  (C-9) 


Equation  (C.9)  is  known  as  the  Radial  Equilibrium  Equation  and  is  a 

linear  differential  equation  with  variable  coefficients. 

If  three  additional  important  assumptions  are  made,  namely  that 

there  is  no  temperature  or  entropy  gradient  at  the  stator  entrance, 
3(RVQ) 

and  that  — =  0  or  free  vortex  flow  exists  ahead  of  the  stator, 

3n  ' 

and  finally  that  the  streamlines  are  cylindrical  surfaces  with  K  =  0, 

3   /   2\ 
then  -r—  (  V   1=0,  which  means  V   is  not  a  function  of  radius.   If 
3n  V  m  J  m 

—  (RV  )  =  0  can  be  maintained  in  the  stator,  the  same  conclusion 

3n    o 


can  be  made  for  V  at  the  stator  discharge, 
m 

The  consequence  of  the  above  series  of  assumptions  results  in  a 
set  of  velocity  profiles  as  shown  in  Fig.  32.   Also  imposed  on  the 
profiles  shown  in  Fig.  32  are  the  effects  of  annulus  friction  and 
secondary  flow. 
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(a)   Axial 

Component  V 


m 


(b)   Tangential 
Component  V 


(c)  Resultant 
Velocity  V 


Figure  32 
Predicted  Velocity  Profiles  at  Stator  Discharge 


The  unequal  pressure  between  the  convex  and  concave  sides  of  a 
blade  cause  trailing  vorticies  to  be  formed  and  are  the  origin  of 
the  secondary  flow.   The  vorticies  cause  a  reduction  in  the  tangen- 
tial component  near  the  hub  and  tip  radii.   It  is  evident  from  the 
above  velocity  profile  in  Fig.  32C  that  the  pressure  distribution 

must  necessarily  be  a  complicated  one  and  not  linear. 

3(RVQ) 


When  the  assumption  that 


nates,  so  that 


3(RV  )      3(RV  tana) 

1.1  3. 


3n 


=  0  is  put  in  cylindrical  coordi- 


9R 


3R 


or 


hence 


V   R  ^tana)   +  y   t£ma 
a     3R        a 


=  0 


=  0 


3 (tana) 

3R 


tana 
R 


(CIO) 
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3(RV  ) 
But  for  — — =  0  the  blades  must  themselves  be  radial,  and  their 

geometry  is  given  by  Vavra  [Ref.  9]  as 


3a     tana  ,J  , ..  * 

M-+TT-  (c.ii) 

Since  equations  (CIO)  and  (C.ll)  are  in  opposition  to  each  other, 

the  blades  cannot  everywhere  be  radial,  and  thus  the  blades  will 

exert  a  force  on  the  flow  in  a  radial  direction,  which  will  produce 

curved  stream  surfaces.   As  a  result,  the  assumption  of  axisymmetric 

3(RV  ) 
stream  surfaces  is  in  error,  and  in  addition  — rr =  0  cannot  be 

OR 

maintained  in  the  stator,  so  free  vortex  flow  is  not  possible.   Thus 

2 
3V         V   3  (RV  ) 

a      »  _u  u    and  this  can  be  solved  if  it  is  assumed  there 

3R         R    3R 

is  constant  total  enthalpy  and  flow  angle  ahead  of  the  stator.   The 

solution  is  given  by  Vavra  [Ref.  8]  which  is: 


where  V   =  axial  velocity  at  the  mean  raduis .   Thus  V  =  f(R)  and 
am  y  a 

the  resulting  velocity  or  pressure  distribution  becomes  extremely 
complicated. 
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